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VECTOR ALGEBRA PACKAGE 


DECUS Program Library Write-up 


DECUS No. 8-254 


ABSTRACT 

This package is designed to operate with the Basic Three Word Floating Point Package 
(Digital 8-5-S) and enables the user to label three dimensional vectors (i.e. three consecutive 
three-word floating point numbers) by a single symbol. The use of this package is similar to 
that of the Floating Point Package in that vector operations can be initiated by a single 
instruction. The operation of addition, subtraction, dot product, cross product and modulus 
of vectors can be simply programmed with this package, e.g. the operation D = (AxJJ)xC 
would be programmed as: 

VGET A 
VCR OS B 
VCROS C 
VPUT D 

Full input and output facilities for vectors are available with this package. 

REQUIREMENTS 

Storage 

This package occupies locations 2-4, 64-74, 5200-5577. The three-word Floating Point 
Package (Basic) occupies locations 5-7, 40-63, 5600-7577. 

Total storage requirements are 2-7, 40-74, 5200-7577. t 

Equipment 

Standard PDP-8, ASR-33 teletype. 

USAGE 

Loading 

This package is loaded with the binary loader (Digital 8-2-U). The binary tape supplied con¬ 
tains only the Vector Algebra Package and hence the three-word Floating Point Package 
(Digital 8-5-S) must also be loaded separately. This package can also be loaded with the EAE 
version of the Floating Point Package, thus reducing the execution times of the vector instruc¬ 
tions. Note that the Vector Algebra Package cannot be loaded with Floating Point Package 
with Extended Functions, since the same area of core is used by each package. 










Calling Sequence 

Location 4 is used as an indirect address for entry to the Vector Algebra Package. Hence 
JMS I 4 will enter the package. After this instruction all subsequent instructions must be of 
the "vector form" (see Discussion) until the exit instruction VEXT is encountered. 

DESCRIPTION 

Discussion 

In the following discussion it is assumed that the reader is familiar with the Floating Point 
Package. The Vector Algebra Package uses locations 64-74 as a vector accumulator, 
similar to the way in which the Floating Point Package uses location 44-46 as a floating 
point accumulator. 


Op. Code 


Mnemonic 


VPUT 


VMOD 


Effect 

Vector Put. 

Store the contents of the vector 
accumulator in the locations specified 
by the effective address. The vector 
accumulator is unchanged. 

Vector Modulus. 

Compute the modulus or length of the 
vector in the vector accumulator and 
store result in the floating point 
accumulator. 


VEXT 


Exit from the Vector Algebra Package. 


Direct or indirect addressing can be used with these instructions. Note that since these 
instructions are not contained in the symbol tables of the PAL or MACRO compilers they must 
be defined by the user's program as follows: 


VADD 

= 

1000 

VSUB 

= 

2000 

VDOT 

= 

3000 

VCR OS 

= 

4000 

VGET 

= 

5000 

VPUT 

= 

6000 

VMOD 

= 

7000 

VEXT 

= 

0000 


Input and Output 

Locations 2 and 3 are used for indirect addressing of the vector input and output routines in 
a manner similar to the way in which locations 5 and 6 are used by the Floating Point Package. 
Hence the instruction JMS I 2 enables the user to input the three components of a vector via 
the teletype into the vector accumulator. 
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Similarly the instruction JMS I 3 enables the user to output the three components of the 
vector contained in the vector accumulator. The three components are output on a single 
line with 5 spaces between each component followed by a carriage-return/line-feed. 

Since the vector input and output routines utilize the input and output routines of the Floating 
Point Package the format of each component is that of a floating point number. Note also 
that the input and output routines use the instructions: 

TSF 

JMP .-1 
TLS 

Hence the teletype flag must be set before input or output. 



Location 


Vector accumulator 

64 

65 

66 

67 

70 < 

71 

72 

73 

74 

r Exponents 
Mantissa > 
Mantissa ) 
Exponents 
Mantissa l 
Mantissa / 
Exponents 
Mantissa > 
a Mantissa ■* 

X - component 

Y - component 

Z - component 

Each component is stored in the normal three-word floating point format. All vectors are 
assumed to be stored in this format where the identifying symbol refers to the location which 
contains the exponent of the X - component. The user is responsible for ensuring that each 
vector is allocated sufficient core storage, i.e. 9 locations. 

Vector Instructions 




Op. Code 

Mnemoni 

c 

Effect 

1 

VADD 


Vector Addition. 

Add vectorially the contents of the 
effective address to the vector accumulator. 

2 

VSUB 


Vector Subtraction. 

Subtract vectorial ly the contents of the 
effective address from the vector accumulator 

3 

VDOT 


Vector Dot Product. 

Form the dot product of the contents of the 
effective address with the vector accumu¬ 
lator. Result in floating point accumulator. 
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Op. Code 
4 


Mnemonic 


VCR OS 


VGET 


Effect 

Vector Cross Product. 

Form the cross product of the vector 
accumulator with the contents of the 
effective address. 

Vector Get. 

Load the vector accumulator with the 
contents of the effective address. 


Example 

As an example of the use of this package consider the calculation of: 

= (AxB)xC 

where A, B, C are vectors to be input via the teletype and D_ is to be output on the teletype. 
/EXAMPLE OF PROGRAM TO COMPUTE D=(AXB)XC 


A=300 


B=311 


C=322 


D=333 


*200 


CLA 


TLS 

/SET FLAG 

JMS 1 2 

/INPUT A 

JMS 1 4 


VPUT A 


VEXT 


JMS 1 2 

/INPUT B 

JMS 1 4 


VPUT B 


VEXT 


JMS 1 2 

/INPUT C 

JMS 1 4 


VPUT C 


VGET A 


VCROS B 

/AXB 

VCR OS C 

/ (AXB)XC 

VPUT D 

/STORE RESULT IN D 

VEXT 


JMS 1 3 

/OUTPUT D 

HLT 
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METHODS 


Discussion 


Let A = ai + a i + ak; B = bi + b j + b k 
— x— yJ- z— — x— yJ- 2 — 


where t_, j_, j< are unit vectors. 

The vector operations are performed as follows: 


Vector Addition 
Vector Subtraction 
Dot Product 


A+B - (a^+b^)j_ + (a^+b^)|_ + (a +b )k 
A-B = (a -b )i + (a -b )j + (a -b )k 

— x X— y y'-L ' 2. z !— 

A.B = a b + ab + a b 
-xx y y z z 


Cross Product 


Vector Modulus 


AxB = (a b -b a )i + (a b -a b )i 
- y z y z — z x x z •*- 


+ (a b -b a )k 
x y x y — 


A I =/a ^ + a 
- 1 v/ X 


2 ^ 2 
+ a 

y z 



5 







/VECTOR PROGRAMMING PACKAGE 


/DEFINITIONS 

SQROOT=2 

FADD=1#0# 

FSUB=^000 

FMPY=3### 

FDIV=4### 

FGET=5### 

FPUT=6^00 

FNOR=7### 

FEXT=#### 


/VECTOR ACCUMULATOR 
VI = 64 /X-CPT. 

V2 = 67 /Y-CPT. 

V3 = 72 /Z-CPT. 



##02 5276 
###3 5314 
0004 5200 
0005 7400 
0006 7200 
0007 5600 

*5200 

5200 0000 VP NT, 

5201 7300 

5202 1600 

5203 2200 

5204 7450 

5205 5600 

5206 3252 

5207 1252 

5210 0254 

5211 7650 

5212 5215 

5213 1 257 

5214 0200 

5215 3261 

5216 126# 

5217 ,0252 

5220 1261 

5221 3261 

5222 1255 

5223 #252 

5224 765# 

5225 523# 


VIN 

VOUT 

VPNT 

74## 

7200 

5600 

0 

CLA CLL 
TAD I VPNT 
ISZ VPNT 
SNA 

JMP I VPNT 
DCA JUMP 
TAD JUMP 
AND PAGENO 
SNA CLA 
JMP .+3 
TAD MASK5 
AND VPNT 
DCA AD DR 
TAD MASK7 
AND JUMP 
TAD ADDR 
DCA ADDR 
TAD INDRCT 
AND JUMP 
SNA CLA 
JMP .+3 


/ADDRESS OF INPUT ROUTINE 
/ADDRESS OF OUTPUT ROUTINE 
/ADDRESS OF PACKAGE 


/GET NEXT INSTRUCTION 


/SET PAGE BIT 
/PAGE ZERO? 
/YES 
/NO 


/GET 7 BIT ADDRESS 
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5226 1661 


TAD 1 ADDR 

5227 3261 


DCA ADDR 

523,0 1261 


TAD ADDR 

5231 3663 


DCA 1 IU1 

5232 1663 


TAD 1 IU1 

5233 1262 


TAD CONF 

5234 3664 


DCA 1 IU2 

5235 1664 


TAD 1 IU2 

5236 1262 


TAD CONF 

5237 3665 


DCA 1 IU3 

5240 1252 


TAD JUMP 

5241 7106 


CLL RTL 

5242 7006 


RTL 

5243 0256 


AND MASK3 

5244 1266 


TAD TABLE 

5245 3253 


DCA JUMP2 

5246 1653 


TAD 1 JUMP2 

5247 3253 


DCA JUMP2 

5250 4653 


JMS 1 JUMP2 

5251 5201 


JMP VPNT+1 

5252 0000 

JUMP, 

0 

5253 0000 

JUMP2, 

0 

5254 0200 

PAGE NO, 

200 

5255 0400 

INDRCT, 

400 

5256 0017 

MASK3, 

17 

5257 7600 

MASK5, 

7600 

5260 0177 

MASK7, 

177 

5261 0000 

ADDR, 

0 

5262 0003 

CONF, 

3 

5263 5564 

IU1, 

U1 

5264 5565 

IU2, 

U2 

5265 5566 

IU3, 

U3 

5266 5255 

TABLE, 

TABLE 

5267 5400 


VAD 

5270 5415 


VSU 

5271 5432 


VDT 

5272 545,0 


VCS 

5273 5521 


VGT 

5274 5533 


VPT 

5275 5545 


VMD 

5276 0000 

VIN, 

0 

5277 4405 


JMS 1 5 

5300 4407 


JMS 1 7 

5301 6064 


FPUT VI 

5302 0000 


FEXT 

5303 4405 


JMS 1 5 

5304 4407 


JMS 1 7 

5305 6067 


FPUT V2 


/INDIRECT ADDRESS OF X-CPT 


/INDIRECT ADDRESS OF Y-CPT 


/INDIRECT ADDRESS OF Z-CPT 


/SELECT OPERATION 


/TABLE FOR INTERPRETING 
/PSEUDO-INSTRUCTIONS 


/ROUTINE TO INPUT A VECTOR 
/INTO VECTOR ACCUMULATOR 

/INPUT X-CPT 


/INPUT Y-CPT 
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53/6 

0000 


FEXT 

53/7 

44/5 


JMS 1 5 

531/ 

44/7 


JMS 1 7 

5311 

6/72 


FPUT V3 

5312 

//// 


FEXT 

5313 

5676 


JMP 1 VIN 

5314 

0000 

VOUT, 

/ 

5315 

7200 


CLA 

5316 

1/55 


TAD 55 

5317 

3361 


DCA SAVE 

53?0 

3/55 


DCA 55 

5321 

44/7 


JMS 1 7 

5322 

5/64 


FGET VI 

5323 

0000 


FEXT 

5324 

44/6 


JMS 1 6 

5325 

4343 


JMS SP 

5326 

44/7 


JMS 1 7 

5327 

5/67 


FGET V2 

533/ 

//// 


FEXT 

5331 

44/6 


JMS 1 6 

5332 

4343 


JMS SP 

5333 

2055 


ISZ 55 

5334 

44/7 


JMS 1 7 

5335 

5/72 


FGET V3 

5336 

0000 


FEXT 

5337 

44/6 


JMS 1 6 

534/ 

1361 

V 

TAD SAVE 

5341 

3/55 

DCA 55 

5342 

5714 


JMP 1 VOUT 

5343 

//// 

SP, 

0 

5344 

1357 


TAD NMB 

5345 

3356 


DCA CTR 

5346 

136/ 


TAD SPACE 

5347 

6/41 


TSF 

535/ 

5347 


JMP .-1 

5351 

6/46 


TLS 

5352 

2356 


ISZ CTR 

5353 

5347 


JMP .-4 

5354 

73/0 


CLA CLL 

5355 

5743 


JMP 1 SP 

5356 

0000 

CTR, 

/ 

5357 

7772 

NMB, 

-6 

536/ 

0240 

SPACE, 

24/ 

5361 

0000 

SAVE, 

/ 


/INPUT Z-CPT 


/ROUTINE TO OUTPUT A VECTOR 
/FROM VECTOR ACCUMULATOR 


/OUTPUT X-CPT 


/OUTPUT Y-CPT 


/OUTPUT Z-CPT 


/ROUTINE TO INSERT SPACES 
/BETWEEN CPTS 
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/VECTOR ROUTINES 




540/ 

//// 

/VECTOR ADD 
*54/0 

VAD, / 

54/1 

44/7 

JMS 1 7 

54/2 

5/64 

FGET VI 

54/3 

1764 

FADD 1 U1 

54/4 

6/64 

FPUT VI 

54/5 

5/67 

FGET V2 

54/6 

1765 

FADD 1 U2 

54/7 

6/67 

FPUT V2 

541/ 

5/72 

FGET V3 

5411 

1766 

FADD 1 U3 

5412 

6/72 

FPUT V3 

5413 

//// 

FEXT 

5414 

56/0 

JMP1VAD 

5415 

//// 

/VECTOR SUBTRACT 

VSU, / 

5416 

44/7 

JMS 1 7 

5417 

5/64 

FGET VI 

54?0 

2764 

FSUB 1 U1 

5421 

6/64 

FPUT VI 

5422 

5/67 

FGETV2 

5423 

2765 

FSUB 1 U2 

5424 

6/67 

FPUT V2 

5425 

5/72 

FGET V3 

5426 

2766 

FSUB 1 U3 

5427 

6/72- 

FPUT V3 

543/ 

//// 

FEXT 

5431 

5615 

JMP 1 VSU 

5432 

//// 

/VECTOR DOT PRODUCT 
VDT, / 

5433 

44/7 

JMS 1 7 

5434 

5/64 

FGET VI 

5435 

3764 

FMPY 1 U1 

5436 

6/64 

FPUT VI 

5437 

5/67 

FGET V2 

544/ 

3765 

FMPY 1 U2 

5441 

6/67 

FPUT V2 

5442 

5/72 

FGET V3 

5443 

3766 

FMPY 1 U3 

5444 

1/67 

FADD V2 

5445 

1/64 

FADD VI 

5446 

//// 

FEXT 

5447 

5632 

JMP 1 VDT 
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/VECTOR CROSS PRODUCT 


545/ 

0000 

VCS, 

/ 

5451 

44/7 


MS 1 7 

5452 

5/72 


FGET V3 

5453 

3765 


FMPY 1 U2 

5454 

63/5 


FPUT TEMPI 

5455 

5/67 


FGET V2 

5456 

3766 


FMPY 1 U3 

5457 

23/5 


FSUB TEMPI 

546/ 

63/5 


FPUT TEMPI 

5461 

5/64 


FGET VI 

5462 

3766 


FMPY 1 U3 

5463 

6311 


FPUT TEMP2 

5464 

5/72 


FGET V3 

5465 

3764 


FMPY 1 U1 

5466 

2311 


FSUB TEMP2 

5467 

6311 


FPUT TEMP2 

547/ 

5/67 


FGET V2 

5471 

3764 


FMPY 1 U1 

5472 

6315 


FPUT TEMP3 

5473 

5/64 


FGET VI 

5474 

3765 


FMPY 1 U2 

5475 

2315 


FSUB TEMP3 

5476 

6/72 


FPUT V3 

5477 

5311 


FGET TEMP2 

55// 

6/67 


FPUT V2 

55/1 

53/5 


FGET TEMPI 

55/2 

6/64 


FPUT VI 

55/3 

//// 


FEXT 

55/4 

565/ 


JMP 1 VCS 

55/5 

//// 

TEMPI, 

/ 

55/6 

//// 


/ 

55/7 

0000 


/ 

551/ 

//// 


/ 

5511 

//// 

TEMP2, 

/ 

5512 

0000 


/ 

5513 

0000 


/ 

5514 

//// 


/ 

5515 

//// 

TEMP3, 

/ 

5516 

000 


/ 

5517 

0000 


/ 

5520 

000 


/ 


/TEMPORARY STORES USED 
/BY CROSS PRODUCT ROUTINE 
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V 


/VECTOR GET 



5521 //// 

VGT, 

/ 


5522 44/7 


JMS 1 7 


5523 5764 


FGET 1 U1 


5524 6/64 


FPUT VI 


5525 5765 


FGET 1 U2 


5526 6/67 


FPUT V2 


5527 5766 


FGET 1 U3 


553/ 6/72 


FPUT V3 


5531 //// 


FEXT 


5532 5721 

/VECTOR PUT 

JMP1VGT 


5533 //// 

vpt. 

/ 


5534 44/7 


JMS 1 7 


5535 5/64 


FGET VI 

% 

5536 6764 


FPUT 1 U1 


5537 5/67 


FGET V2 


554/ 6765 


FPUT 1 U2 


5541 5/72 


FGET V3 


5542 6766 


FPUT 1 U3 


5543 //// 


FEXT 


5544 5733 


JMP 1 VPT 



/VECTOR LENGTH 


5545 //// 

VMD, 

/ 


5546 44/7 


JMS 1 7 


5547 5/64 


FGET VI 


555/ 3/64 


FMPY VI 


5551 6/64 


FPUT VI 


5552 5/67 


FGET V2 


5553 3/67 


FMPY V2 


5554 6/67 


FPUT V2 


5555 5/72 


FGET V3 

# 

5556 3/72 


FMPY V3 


5557 1/67 


FADD V2 


556/ 1/64 


FAD VI 


5561 ///2 


SQROOT 


5562 //// 


FEXT 


5563 5745 


JMP 1 VMD 


5564 //// 

ui. 

/ 


5565 //// 

U2, 

/ 


5566 //// 

U3, 

/ 


/USED AS INDIRECT ADDRESS 
/OF OPERAND 
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ADDR 

5261 

CONF 

5262 

CTR 

5356 

FADD 

1000 

FDIV 

4000 

FEXT 

0000 

FGET 

5000 

FMPY 

3000 

FNOR 

7000 

FPUT 

6000 

FSUB 

2000 

INDRCT 

5255 

IU1 

5263 

IU2 

5264 

IU3 

5265 

JUMP 

5252 

JUMP2 

5253 

MASK3 

5256 

MASK5 

5257 

MASK7 

5260 

NMB 

5357 

PAGE NO 

5254 

SAVE 

5361 

SP 

5343 

SPACE 

5360 

SQROOT 

0002 

TABLE 

5266 

TEMPI 

5505 

TEMP2 

5511 

TEMP3 

5515 

U1 

5564 

U2 

5565 

U3 

5566 

VAD 

5400 

VCS 

5450 

VDT 

5432 

VGT 

5521 

VIN 

5276 

VMD 

5545 

VOUT 

5314 

VPNT 

5200 

VPT 

5533 

VSU 

5415 

VI 

0064 

V2 

0067 

V3 

0072 




# 
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